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Abstract. This work deals with first hitting time densities of Ro 
processes whose local drift can be modeled in terms of a solution to 
Burgers equation. In particular, we derive the densities of the first 
time that these processes reach a moving boundary. We distinguish 
two cases: (a) the case in which the process has unbounded state 
space before absorption, and (b) the case in which the process has 
bounded state space before absorption. The reason as to why this 
distinction has to be made will be clarified. 

Next, we classify processes whose local drift can be expressed as 
a linear combination to solutions of Burgers equation. For example 
the local drift of a Bessel process of order 5 can be modeled as the 
sum of two solutions to Burgers equation and thus will be classified 
as of class B 2 . Alternatively, the Bessel process of order 3 has a 
local drift that can be modeled as a solution to Burgers equation 
and thus will be classified as of class B 1 . Examples of diffusions 
within class B 1 , and hence those to which the results described 
within apply, are: Brownian motion with linear drfit, the 3D Bessel 
process, the 3D Bessel bridge, and the Brownian bridge. 



I. Introduction 

This work deals with processes whose local drift is modeled in terms 
of a solution to Burgers equation 

-/Hit, x) = ^fi xx (t, x) + /i(t, x) ■ fx x (t, x), (t, x) eR + x R. 
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Some of these processes appear prominently in stochastic analysis. Its 
applications in finance, as well as its appearance in statistical prob- 
lems is also noteworthy Among examples of processes with this prop- 
erty we can mention: Brownian motion, Brownian motion with lin- 
ear drift [for financial applications of this process see Section 5.5.8 in 
Karatzas and Shreve (1991)], the 3D Bessel process, the 3D Bessel 
bridge [some recent financial applications of this process can be found 
in Davis and Pistorius (2010), it is also related to the Kolmogorov- 
Smirnov and Cramer von Mises tests, see Gikhman (1957) and Kiefer 
(1959)], and the Brownian bridge [for applications of this process in 
tests of Kolmogorov-Smirnov type see Kolmogorov (1933), Smirnov 
(1948), and Andel (1967)] . 

In fact any Bessel process of odd order greater than or equal to 3 has 
a local drift which can be modeled as linear combinations of solutions 
to Burgers equation. For instance, the Bessel process of order 5 has 
a local drift which is a sum of two solutions of Burgers equation. In 
general a Bessel process of order 2k + 1, k > 1, has a local drift which 
can be written as a sum of k solutions to Burgers equation. This leads 
us to classify SDEs to be of class B k if its local drift is a sum of k 
solutions to Burgers equation. Regarding the related topic of hitting 
times of Bessel processes there are a number of recent papers: e.g. 
Salminen and Yor (2011) or Alili and Patie (2010). 

In this paper we derive the densities of the first time that processes 
of class B l reach moving boundaries. Furthermore we assume that 
the moving boundaries are real valued and twice continuously differen- 
tiable. To this end we distinguish two cases: (a) the case in which the 
process has unbounded state space before absorption, and (b) the case 
in which the process has bounded state space before absorption. An 
example of the second case is the density of the first time that a 3-D 
Bessel bridge started at y > 0, and absorbed at zero at time s, hits a 
fixed level a, where y < a. That is, the 3-D Bessel bridge lives on (0, a) 
before being absorbed at either level a, or at level at time s, see for 
instance Hernandez-del- Valle (2012). 

The paper is organized as follows in Section [2] an /i-transform and 
general notation are introduced. Next, in Section |3j the unbounded 
state space case is discussed. Section |4j deals with problems where the 
process lives on a bounded interval. Next, in Section |5j we give a more 
detailed description of class B 1 and its relationship with class £> 2 .We 
conclude in Section |6] with some final comments and remarks. 
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2. Preliminaries 

Remark on notation 2.1. As in the analysis of diffusion processes, 
PDEs with derivatives with respect to (t, x) are called backward equa- 
tions, whereas PDEs with derivatives in (s, y) are called forward equa- 
tions. 

Furthermore, through out this work (i) B = {B t , J~t}t>o stands for 
one- dimensional standard Brownian motion, (ii) For a given function, 
say w, partial differentiation with respect to a given variable, say x, 
will be denoted as w x . 

In the following sections our main tools will be: (1) An /i-transform, 
see Theorem 2.2, and (2) Ito's lemma, see Lemma 2.3[ [Regarding 
/t-transforms the reader may consult Doob (1984) or Pinsky (1995).] 

Theorem 2.2. Let h be of class C 1,2 (R+ x R) as well as a solution to 
the backward heat equation 

(1) -ht= 7^h xx . 

Furthermore, consider processes X , and Y which respectively satisfy 
(at least in the weak sense), the following equations (each under their 
corresponding measures P, and Q), 

(2) (P) dx t = h ^4rdt + dB t 

n(t, A t ) 

(3) (Q) dY t = dB t . 

Moreover, suppose that f is real valued and integrable. Then the fol- 
lowing identity holds 



(4) KMX r )\ = Eg, 



n{t, x) 



Proof. For the proof, see Theorem 2.10 in Hernandez-del- Valle (2011). 

□ 

We recall Ito's lemma in the case in which the drift is a deterministic 
function of time. 

Lemma 2.3. Let /(•) be a real-valued differentiable function, h a so- 
lution of the one- dimensional backward heat equation, and processes Z 
and S have the following dynamics 

dZ t = f'(t)dt + dB t 
dS t = -f(t)S t dB t 
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for < t < oo, under some meaure Q. Then 

S.-h(-,Z.) 

is a Q-martingale. 
Proof. From Ito's lemma 

dh(t,Z t ) = ht(t,Z t )dt + h z (t,Z t )f'(t)dt 

+h z (t,Z t )dY t + -h zz (t,Z t )dt 
= h z (t,Z t )f(t)dt + h g (t,Z t )dY t . 

Hence, 

d[S t -h(t,Z t )\ = f(t,Z t )dS t + S t dh(t,Z t ) 

+dh{t, Z t ) ■ dS t 
= h{t,Z t )[-f'{t)S t dY t ] 

+S t h g (t, Z t )f'(t)dt + S t h z (t, Z t )dY t 

-h z (t,Z t )S t f(t)dt 
= [S t h z (t,Z t )-f'(t)X t h(t,Z t )]dY t . 

3. Unbounded state space 



□ 



Remark 3.1. Henceforth let h be a solution to the one- dimensional 
backward heat equation Let process X have the following dynamics 



,„ h x (t,X t ) f R or 

(5) dX t = — — dt + dB t , X = y e < R+ 

for < t < oo. And assume the drift h x /h satisfies the Ito conditions. 

In this section, we specialize to the case in which the process X, with 
dynamics as in ([5]), has unbounded state space before absorption. In 
particular we find the density of the first time that X hits a real- valued 
and twice continuously differentiable function /. To this end let us first 
define the following stopping times. 

Definition 3.2. Given the constant a G K and the real-valued, twice 
continuously differentiable function /(•) — which we refer to as a "mo- 
ving boundary" — we define the following stopping times 

(6) T := inf jt > 0|X t = a + jf /'(u)duj y^a 
T B := inf jt > 0\B t = a + J /'(u)duj . 
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Furthermore, letpf(-) be the density of T B , that is, the density of the 
first time that a one- dimensional Wiener process reaches a determin- 
istic moving boundary f . For detailed historical and technical account 
(making use of integral equations) of this problem see Peskir (2001). 

Next, we present the main result of this section. 

Theorem 3.3. Suppose that X has dynamics as in T is as in 
and a ^ y. Then 



F y (T E du) 



h (u,a + J Q U f'(v)dv) 
h(0,y) 



pf(u)d 



u. 



for u > . 

Proof. For T as in (JSj) , it follows from Theorem 2.2 that 
P y (T<t) = E y [l (T<t) ] 
h(t,Y t ) 



E' 



h{0,Y o 



-I 



(T<t) 



From Girsanov's theorem, and given that Y is a Q- Wiener process we 
have 

- /« f'(u)dY u -i [:(f(u)^du h ( t ^ Yt + Jq f\u)du) ^ 

h(0,y) 



E 



-I 



(T<t) 



Finally, from Lemma 2.3 and the optional sampling theorem 

- rJ fw.-h tfifM)** h ( T > a + Jo f'( u ) du ) * A 

h(0,y) (T<t) 



E l 



h[u,a + J Q U f'{y)dv) 
h(Q,y) 



pf(u)d 



u. 



□ 



Examples of processes which satisfy equation ^ and that have un- 
bounded domain before hitting the boundary / are: 

Examples 3.4. (i) Standard Brownian motion, where h(t, x) = c. 
(ii) Brownian motion with linear drift, where 



h(t, x) = e" 
(iii) Brownian bridge, where 
1 



x 



h(t, x) 



_g 2(s— *) 



(iv) 3D Bessel process, where h(t,x) = x. 



(t, x) e [0, s] x 



6 GERARD O HERNANDEZ-DEL-VALLE 

(v) 3D Bessel bridge, where 

hit.x) = X - e ~^), (t, x) G [0, s) x ]R + . 

Example 3.5. In Figures^ and\^ the theoretical density and distribu- 
tion of the first time that a Brownian bridge started at y = 1, absorbed 
at time s = 3 at level c = hits the linear barrier 

f(t) = 2-t, t>0 

is compared with n = 5500 simulations. See Durbin and Williams 
(1992). 

Example 3.6. In Figures^ and\^ the theoretical densities and distri- 
butions of the first time that a 3-D Bessel process started at y = 3 and 
absorbed at time s = 4, reaches level a — 1. For a general overview 
of the 3D Bessel bridge, see Revuz and Yor (2005). The case in which 
a level is reached from below is in general studied in Pitman and Yor 
(1999) or Hernandez- del- Valle (2012). 

4. Bounded domain 

In this section we will also also use the following: 

Given that B is a one-dimensional Wiener process started at y and 

T := inf{t >0\B t = 0} 
T := inf {t > 0\B t = a}, < y < a. 



We recall 



P„(T A T G dt) := -== V 



(2na + y) exp 
+ {2na + a — y) exp 



(2na + yf 



2t 

(2na + a — y) 2 



F tv (T G s,T > s) := V 



n=—oo 



2t 

{2na + a — y) 
(2na + a — y) 2 



dt, 



2(s - t) 



x exp 

P„(T G dt) := , V exp {- — 1 dt. 
yK J V2^ I 2t) 

See for instance Chapter 2, Section 8 in Karatzas and Shreve (1991). 
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In this section, we specialize to the case in which the process X, with 
dynamics as in (|5]), has bounded state space before absorption. The 
main result of this section is the following. 



Theorem 4.1. Given that X has dynamics as in pD and T is as 

We have for < t < s and < y < a that 

F y (T G du) = ; ' , ; [P y (T A T G du) - ¥ y (T G du)} , 



in 



h(0,y) 



where < u < s. 



Proof. We follow the proof of Theorem 3.3 However, we must take 
into account the fact that the Q- Wiener process Y is absorbed at zero. 
Given that T is as in ^ 

Py(T<t) = E y [l (T<t) I (2b>t) ] 

h{t,Y t ) 



E 



E 



h(o,Y y (T<t > To>t \ 

h( T , a) tt 

H^Yo) kT<t ' To>t \ 



where the last line follows from the optional sampling theorem. Finally 
recall the identity 

F y (T<t,T >t) = F y (T >t)-F y (T>t,T >t) 
= Fy(T >t)-F y (TAT >t) 
= ¥ y (T AT <t)-¥ y (T <t). 

□ 

Examples 4.2. Some examples are: (i) 3D Bessel process (reaching 
a fixed level from below) (ii) 3D Bessel bridge (reaching a fixed level 
from below). Regarding the first hitting probabilities of general Bessel 
processes see Wendel (1980), or Betz and Gzyl (1994a, 1994b). 

Example 4.3. In Figure\5\ the theoretical distribution of the first time 
that a 3D Bessel process reaches a = 1.5 from below, is compared with 
a simulation n = 5500 (see hard line). 

Example 4.4. If we set 

h a (s-t,x) :=F tjX (T es,T > s), 



GERARDO HERNANDEZ-DEL- VALLE 



and define a process Y to be as in 

h a y (s-t,Y t ) 



dY t = ^dt + dB t , 0<t<s 

h a (s-t,Y t ) 



One can show that this process has state space (0, a) for t G [0, s) [see 
Hernandez- del- Valle (2011)]. In Figure^the density and distribution 
of the Wiener process, started at x = 1/2, absorbed at zero and that 
reaches level a = 2 for the first time at s = 2 is plotted at t = 1 and 
£ = 7/8. 

5. Heat polynomials and Burgers equation 

In this section we provide a characterization of the results described 
in Sections |3] and HI 

To this end let us first introduce the following classification of SDEs. 

Definition 5.1. We will say that process X , which satisfies the follow- 
ing equation 

dX t = fi(t, X t )dt + dB t , 
is of class B n , n = 1,2, . . . , if its local drift \x can be expressed as 

K ' v ' ^ h>(t,x) 

Where each W is a solution to the backward heat equation |7]). 

Making use of this classification it follows that: 

Remark 5.2. Processes X\, X2, and X3 which respectively satisfy the 
following equations 

[F 1 ) dXt(t) = jfjfidt + dB t 

;P 2 ) dX 2 (t) = ~^dt + dB t 



(P 3 ) dX 3 {t) 



dt + dB t 



X 3 (t) s-t 

are of class B . That is, the 3-D Bessel process X\, the Brownian bridge 
X2, and the 3-D Bessel bridge have a local drift which is a solution to 
Burgers equation. This statement is verified by using the following 
solutions to the heat equation correspondingly 



(9) 



k(t,x)=x, g {t,x) = -j^= T) e W {-^) 

^ = 7^ exp {-2(fe)}' 
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together with the Cole-Hopf transform which relates Burgers equation 
with the heat equation. 

Examples of processes which are not B l are the following. 



Example 5.3. The Bessel process of odd order 2n + 1, n = 1, 2, . 
of class B n . 

Recall that the Bessel process of order m G N is the solution to 

1 



■is 



2X t 



'-dt + dB t 



if m = In + 1 
dX t 



(2n + l) 



1 



2X t 
-dt + dB t 



dt + dB t 



n 



k x (t,X t ) 
k(t,X t ) 



+ ... _|_ 



k x (t,X t 
k(t,X t ) 



dt + dB t 



where k is as in |$|). 

However there exists a least one important process which is both B 1 
and B 2 . 

Proposition 5.4. The 3-D Bessel bridge process X 3 , which solves 
is B 1 and B 2 . 



Proof, it follows by verifying that for k, g, and h as in @ the following 
identity holds 



h 



9x 



□ 



In turn, this feature of process X 3 [i.e. a process that is both B 1 
and B 2 ], leads to some interesting properties. [See Hernandez-del- Valle 
(2011).] The next natural question to ask is if there are more processes 
with such property. 

In this direction we will show that process X3 is the only element of 
both B l and B 2 . In the case in which the generating functions w, of 
class B 1 , are heat polynomials. [That is if w is a solution to the heat 
equation and w x /w models the local drift of a process X] 

Definition 5.5. (Heat polynomials) [Widder and Rosenbloom (1959)]. 
A heat polynomial Vj (x, t) of degree j is defined as the coefficient of 
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z n /n\ in the power series expansion 

oo 

(10) e xz+\zH = "y v f xt y 



n=0 







An associated function w n (x,t) is defined as 

(11) w n (x,t) = g(t,x)v n (x,-t)(t/2)-'\ 

where g, as in is the fundamental solution to the heat equation. 



Remark 5.6. Observe that if processes X\, X 2 , and X 3 are as in fSp 
then their local drifts can be described in terms of heat polynomials, 
Definition 5.5. In the case of Xi its corresponding polynomial is V\. 
Alternatively for X 2 and X 3 their corresponding polynomials are W\ 
and ui2 respectively. 

The main result of this section is the following. 

Theorem 5.7. If process X 6 B l . [That is X is a solution to 

h x {t,X t ) 



-dt + dB t 



h(t,X t ) 

and h solves the backward heat equation —h t = \h xx .] And h is either 
a heat |7o[ ) or derived heat (11) polynomial. Then the only process 
which is also B 2 is the 3-D Bessel bridge X3, which has dynamics as 
in QF 3 ). 



Proof. Given that v n and w n are as in (10) and (11) respectively; and 
letting w' stand for differentiation with respect to the first variable we 
have that 

w'Jxjt) x v[(x,t) I 

— — and 
t 



Wo(x,t) t Vi(x,t) X 

Next, since [see p. 225 in Rosenbloom and Widder (1959)] 



<-i(M) = --w n (x,t) 



it follows that 



w' n (x,t) x 2nw„_ 1 (x,i) 

w n (x,t) t t w n (x,t) 



(12) 



x n 
~1 ~ I 

w (x,t) 



^(x,t) 



W n 

n 
' 1 



w n _i(x,t) 
w' n _ x (x,t) 
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Alternatively from (11) 



Wn-i _ g{x,t)v n -i(x,-t)(t/2) 



-n+l 



W 



n-1 



-n+l 



-\w n {x,t) 
g(x,t)v n _i(x, -t)(t/2) 
-±g(x,t)v n (x,-t)(t/2)-» 
tv n -i(x, -t) 



v n (x, -t) 



This implies, from (12), that 



(13) 



However, since 



w' n (x,t) 


w' 




t) 


w n (x,t) 


w 




t) 




w' 




t) 




Wo 




t) 




w' 




t) 




Wo 




t) 



t 

n 
1 

+ n 



w ra _i(g,t) 

K-i( x >t) 

tv n _i{x, -t) 
v n (x, -t) 

V n -l(x, -t) 

v n (x,-t) 



v' n (x,t) = nv n -i(x,t) 

[see equation (1.9) in Widder and Rosenbloom (1959)] we have, from 
@, that 

w' n (x,t) w' (x,t) nvn-xjx, -t) 
w n (x,t) w (x,t) v n (x,-t) 

w (x,t) v n (x,-t)' 

In general, since v n (x,t) is a solution to the backward heat equation, 
then v n (x, —t) is a solution to the forward equation. This is true as 
long as n > 1. However if n = 1, vi(x,—t) is also a solution to the 
backward equation becuase it does not depend on t. In this case we 
have that 

w[(x,t) 



Wi_(x,t) 



x 1 
t x 



as claimed. 



□ 



6. Concluding remarks 

In this work we study processes X, which have local drift modeled in 
terms of solutions to Burgers equation. In particular, we find the den- 
sity of the first time that such processes hit a moving boundary. Next, 
we propose a classification of SDEs in terms of solutions of Burgers 
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equation. We say that process X is of class B^ if its local drift can be 
expressed as a sum of j solutions to Burgers equation. 

We note that the 3-D Bessel process, the 3-D Bessel bridge, and 
the Brownian bridge are all B 1 . However we show that the 3-D Bessel 
bridge. Furthermore, we show, that it is the only process which satisfies 
this property if the solutions of Burgers equation is constructed by use 
of heat polynomials. A more detailed study of this classification is work 
in progress. 
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Histogram of 5500 sim. 



Theoretical density 





Figure 1. (Example 3.5). The graph is plotted in R. 
The upper left graph is the histogram of the (simulated) 
first time that a Brownian bridge started at y — 1, and 
absorbed at c = at time s = 3, reaches the linear 
boundary f(t) = 2 — t. In the upper right frame we 
have its theoretical density. In the lower left we have the 
simulated distribution, and finally on its right we have 
its theoretical counterpart. 
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Theoretical vs. Sim. 
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Figure 2. (Example |3.5|). The graph is plotted in R. 
The dotted line is the theoretical probability. The hard 
line is a simulation with n = 5500. 
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Figure 3. (Example 3.6). The graph is plotted in R. 
The upper left graph is the histogram of the (simulated) 
first time that a 3-D Bessel bridge started at y — 3, and 
absorbed at s — 4, reaches level a = 1. In the upper 
right frame we have the its theoretical density. In the 
lower left we have the simulated distribution, and finally 
on its right we have its theoretical counterpart. 
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Theoretical vs. Sim. 
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Figure 4. (Example 3.6). The graph is plotted in R. 
The dotted line is the theoretical probability. The hard 
line is a simulation with n = 4500. 
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Figure 5. (Example 3.6). The graph is plotted in R. 
The dotted line is the theoretical probability with a = 
1.5. The hard line is a simulation with n = 5500. 
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Figure 6. (Example 4.3). The graphs are plotted in R. 



The upper left graph is the density of a Wiener process 
(started at x = 1/2), absorbed at zero and that reaches 
a = 2 for the first time at s = 2; evaluated at t — 1. 
Th upper right graph is the corresponding distribution 
at t — 1. In the lower left figure we have the density 
at time t = 7/4. Finally, the lower right panel is its 
corresponding distribution. 



